Ab initio modeling of Bose-Einstein condensation in Pb 2 V 3 Og 
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We apply density functional theory band structure calculations and quantum Monte Carlo sim- 
ulations to investigate the Bose-Einstein condensation in the spin-i quantum magnet Pb2V30g. 
In contrast to previous conjectures on the one-dimensional nature of this compound, we present a 
quasi-two-dimensional model of spin dimers with ferromagnetic and antiferromagnetic interdimer 
couplings. Our model is well justified microscopically, and provides a consistent description of the 
experimental data on the magnetic susceptibility, high-field magnetization, and field vs. tempera- 
ture phase diagram. The Bose-Einstein condensation in the quasi-two-dimensional spin system of 
Pb2V30g is largely governed by intralayer interactions, whereas weak interlayer couplings have a 
moderate effect on the ordering temperature. The proposed computational approach is an efficient 
tool to analyze and predict high-field properties of quantum magnets. 

PACS numbers: 75.30.Et, 75.10.Jm, 71.20.Ps, 75.50.Ee 



Bose-Einstein condensation (BEC) is one of the fun- 
damental phenomena in physics. Recent activity in the 
field of quantum magnetism opened new prospects for ex- 
tensive experimental investigation of this phenomenon*! 
The bosonic nature of magnons and their control by 
the external magnetic field give rise to simple realiza- 
tion of the BEC, whereas a variety of spin lattices 
available in transition-metal compounds lead to different 
regimes of bosonic interactions. However, the detailed 
understanding of the underlying magnetic couplings re- 
mains an essential prerequisite for the correct interpre- 
tation of the observed high-field physics. For example, 
the proposed two-dimensional (2D) regime of BEC in 
BaCuSi206 (Ref. [2|) was based on an oversimplified spin 
model and later challenged by direct observations of in- 
equivalent spin dimers in this compoundi" A subsequent 
theoretical study reconciled all the experimental observa- 
tions within an extended three-dimensional (3D) model 
comprising two inequivalent sublattices with different bo- 
son densities^ 

To derive a complete spin model, inelastic neutron 
scattering studies are usually required. A viable alter- 
native is given by density functional theory (DFT) band 
structure calculations that are able to evaluate individual 
exchange couplings in specific compounds. Here, we show 
that DFT calculations combined with efficient numerical 
techniques lead to a quantitative description of the BEC 
in Pb2V30g. Based on this approach, the role of individ- 
ual exchange couplings can be explored, and the regime 
of the BEC in a quantum magnet can be predicted. 

Pb2V30g is a model spin-i material showing 
magnetic-field-induced long-range ordering (LRO) inter- 
preted as the BEC of magnons. 6 ' - Despite a number of 
experimental studies available, the microscopic magnetic 
model of this compound remains controversial. The crys- 
tal structure comprises chains of corner-sharing V +4 06 
octahedra ( "structural chains" ) which are joined into lay- 
ers by non-magnetic V +5 04 tetrahedra (Fig. [lj. The 
magnetic behavior is roughly captured by the model of 
alternating spin-i chains. Refs. @ and suggest that 
these spin chains coincide with the structural chains. By 




FIG. 1. (Color online) Crystal structure (top) and spin model 
(bottom) of Pb2V30g. Circles denote positions of V +4 and 
sites of the spin lattice. The shading shows a single spin 
dimer. The right panel depicts a single chain of corner-sharing 
V06 octahedra and respective Wannier functions with d xy 
orbital character. Thick orange lines show the short V-O 
bonds. 



contrast, Mentre et al& found alternating spin chains per- 
pendicular to the structural chains. The latter study also 
proposed ferromagnetic (FM) couplings between the spin 
chains, although no detailed comparison to the experi- 
mental data was given. 

A microscopic investigation of the uniform-spin-chain 
compound Sr 2 V30g with a similar crystal structure^ sug- 
gests that leading antiferromagnetic (AFM) couplings 
run perpendicular to the structural chains. Taking into 
account the lower symmetry (monoclinic in Sr 2 V30g and 
triclinic in Pb2VsOg), one would expect AFM alternat- 
ing spin chains perpendicular to the structural chains and 
thus support the model by Mentre et al& However, the 
couplings between the spin chains are also important. 
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Below, we show that PD2V3O9 should be considered as a 
quasi-2D system in contrast to the quasi-one-dimensional 
(ID) spin lattice in Sr 2 V30g. 

To elucidate the microscopic magnetic model of 
Pb^VaOg, we performed scalar- relativistic DFT calcu- 
lations using the FPLO code 1 ^ and the local density ap- 
proximation (LDA) with the exchange-correlation poten- 
tial by Perdew and Wangj 1 ^ The V +4 -related states were 
further introduced into a multi-orbital Hubbard model 
that allowed to treat strong correlation effects in the V 
3d shell and to evaluate the exchange couplings. The re- 
sults were cross-checked by a local spin-density approx- 
imation (LSDA)+i7 method with an around-mean-field 
double-counting correction (DCC) that accounts for elec- 
tronic correlations in a mean-field fashion. LSDA+C7 to- 
tal energies for collinear spin configurations were mapped 
onto the classical Heisenberg model to obtain individ- 
ual exchange couplings J;. The fc-mesh comprised 512 
points for the LDA calculation and 50 — 100 points for 
the LSDA+J7 supercell calculations. 

The LDA band structure of Pb 2 V 3 9 (Fig. 02) re- 
sembles that of A2V3O9 with A = Sr and Ba£ The 
filled valence bands between —7 eV and —2 eV are 
formed by O 2p states, whereas the bands at the Fermi 
level originate from 3d states of octahedrally coordinated 
vanadium (i.e., V +4 ). Owing to the complete charge 
orderingf&i^ the 3d states of tetrahedrally coordinated 
vanadium atoms (i.e., V +5 ) appear above 1 eV only. Lead 
orbitals give rise to narrow bands around —9 eV (6s) and 
to a pronounced contribution above 2 eV (6p), similar to 
other Pb +2 compounds^ 

V +4 bands show a characteristic crystal-field splitting 
with t 2g states lying below 1 eV and e g states spanning 
from 1 eV to 4 eV. To extract the relevant microscopic 
information, we fit the t 2g bands with a three-orbital 
tight-binding model,™ based on Wannier functions with 
proper orbital characters* 1 ^ The fit reveals a lower en- 
ergy of the d xy orbital (0.04 eV) and higher energies 
of the d yz and d xz orbitals (0.37 and 0.50 eV, respec- 
tively) where z denotes the direction of the shortest V- 
O bond that aligns with the structural chains (Fig. [1} . 
Strong electronic correlations typical of V +4 compounds 
stabilize the half-filled d xy orbital in the Mott-insulating 
state (e.g., via the LSDA+J7 calculation). Therefore, the 
exchange couplings can be estimated using the Kugel- 
Khomskii model; 16 ' 17 

At 2 At 2 1 a 

j _ ^ L xy asy— Va^en 

' " _ a kt, xz (tfcff + A <*)(UeS + A a - J off ) ' U 

where t xy and t xy ^ a are the hoppings between the xy 
states and from the xy (half-filled) to a (empty) states, 
U e s and J e g are the effective on-site Coulomb repulsion 
and Hund's coupling in V 3d bands, respectively, and A Q 
is the crystal- field splitting between the xy and a states. 
The first term in Eq. (TIJ corresponds to AFM couplings 
due to the hoppings between the half-filled xy states. The 
second term is the FM coupling caused by the hoppings 
to the empty states and by the Hund's coupling. 



TABLE I. Exchange integrals (in K) calculated using the 
model approach [Eq. Q] and LSDA+[/ supercell approach, 
see text for details. 
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Using the typical values of U e g = A eV and J c ff = 
1 eVji^ii 7 - we find two AFM and two FM couplings in 
Pb 2 V 3 09 (TableHJ. The stronger AFM couplings J x and 
J[ run between the structural chains via VO4 tetrahcdra, 
whereas the weaker FM couplings J 2 and J 2 run along the 
structural chains (Fig.[T]). This counter- intuitive scenario 
is readily explained by the orbital state of vanadium. In 
terms of the crystal-field theory, the energy preference of 
the xy orbital is caused by the short V-0 bond which is 
directed along the structural chain. The half-filled (xy) 
orbital is therefore located in the plane perpendicular to 
the structural chain (Fig. [I]), and favors the long-range 
couplings Ji and J[. On the other hand, it precludes 
the V-O-V superexchange, despite the V-O-V angles 
are about 145°, i.e., far from 90° (see also Refs. H and 
[Lsh . The resulting couplings along the structural chains 
( J2, J 2 ) are weak and FM due to the short V-V distances. 
A similar scenario has been established for (VO)2P207, 
S^VsOg, and other V +4 compounds ^ 18 ' 19 

Our model basically follows the earlier proposal by 
Mentre et alM and differs from Refs. [y and [j], which 
assume AFM couplings along the structural chains. In 
contrast to the result of Mentre et al. (| J 2 \ and \J' 2 \ are 
comparable to or even larger than J\ and J[)r we find 
rather weak FM couplings that are in agreement with the 
experiment (see below) and do not inhibit the formation 
of the spin gap in Pb 2 V30g. 

The model analysis using Eq. (JIJ is confirmed by 
LSDA+C/ calculations. Employing the on-site Coulomb 
repulsion and exchanged parameters Ud — 3 eV and 
Jd = 1 eV, respectively*^ we find good agreement with 
the results of the model analysis (Tabled]). The estimated 
band gap of 1.55 eV conforms to the reported black 
color of Pb2V309f 1 ^ We also used the fully- localized-limit 
DCC and the generalized gradient approximation for the 
exchange-correlation potential, but these changes in the 
computational procedure did not modify the results qual- 
itatively. 

The resulting spin model of Pl^VsOo is shown in the 
bottom part of Fig. [TJ Although it can be viewed as 
AFM alternating spin chains (along a') that are cou- 
pled ferromagnetically (along 6'), the numerical study 
and the comparison to the experiment emphasize the 
quasi-2D nature of the material (see below). Regard- 
ing the interlayer couplings, we note that each V +4 site 
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FIG. 2. (Color online) Top: LDA density of states for 
PD2V3O9. The Fermi level is at zero energy. Bottom: LDA 
bands at the Fermi level (thin light lines) and the fit with 
the three-orbital tight-binding model (thick dark lines). The 
notation of fc-points is T(0, 0,0), X(0.5,0,0), M(0.5, 0.5, 0), 
Y(0, 0.5, 0), and Z(0, 0, 0.5) in units of the respective recipro- 
cal lattice parameters. 
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Field (T) 

FIG. 3. (Color online) Magnetic susceptibility (top) and high- 
field magnetization (bottom) of Pb2Vs09 (circles) fitted with 
the ID alternating-chain model (dashed line) and the 2D 
Ji — J[ — 3% model (solid line, see also Fig. [TJ. The mag- 
netization curve was measured at 1.5 K, whereas the model 
curves are simulated at T/Ji = 0.05. Experimental data are 
taken from Ref. 0. 



has eight neighbors in the adjacent layers (V-V dis- 
tances of 8.7 — 9.2 A), yet most of these contacts are 
inactive (|Jj| < 0.15 K). The leading interlayer coupling 
J± ~ 0.6 K provides two bonds at each lattice site and 
corresponds to the V-V distance of 9.01 A. Thus, the in- 
terlayer couplings are uniform and unfrustrated. In the 
following, we apply this model to interpret the magnetic 
behavior of Pb2V30g. 

Thermodynamic properties and magnetic ordering 
temperatures were simulated using directed loop algo- 
rithm in the stochastic series expansion representation, 22 
as implemented in the ALPS packaged We used 2D and 
3D finite lattices with periodic boundary conditions. The 
lattice size was set to L x L (2D model) and L x L x L/2 
orLxLxL/4 (3D model) with L < 32 (here, L is the 
number of unit cells; each cell comprises four magnetic 
atoms) 1^ 

Figure [3] shows magnetic susceptibility (x) and high- 
field magnetization (M) of Pb2V3 0g along with the sim- 
ulations for the purely one-dimensional (ID) alternat- 
ing spin chain J\ — J[ model (dashed line: J\ ~ 31 K, 
J[/Ji ~ 0.64, g ~ 1.99). This model underestimates x at 
low temperatures and overestimates the spin gap. More- 
over, the shape of the simulated magnetization curve is 
apparently different from the experimental result while 
the fitted g value exceeds the actual g ~ 1.93 from elec- 
tron spin resonance (ESR). 2 » Thus, the alternating-chain 



model is insufficient to understand the magnetic behavior 
of Pb 2 V 3 9 . 

To improve the ID model, we consider interchain cou- 
plings J2 and J 2- For the moment, we assume a sin- 
gle interchain coupling J2 = J 2 = J 2- The effect of in- 
equivalent J2 and J'2 will be discussed later. Since the 
ID model yields an accurate prediction of the satura- 
tion field, additional (i.e., interchain) couplings should 
be FM. Indeed, we find a perfect fit of the experimental 
data with Ji ~ 31 K, J'J J x = 0.65, J2/J1 = -0.25, and 
g ~ 1.94 (compare to 1.93 from ESR)^ 

The fitted exchange parameters are in remarkable 
agreement with our DFT results. The leading AFM cou- 
pling is Ji = 31 K (37 K in LSDA+t/, see Tabled, 
whereas the interchain couplings are indeed FM. The in- 
terchain coupling J2/J1 = —0.25 is sufficient to reduce 
the spin gap of Pb2V30g and smoothen the magnetiza- 
tion curve with respect to the simulation result for the ID 
model. A stronger interchain coupling (as proposed by 
Mentre el a/4) will close the spin gap, thus contradicting 
the experimental magnetic behavior of Pb2V30g. 

The key feature of Pb2V30g is the field-induced LRO 
(i.e., BEC) between the critical field H c and the satu- 
ration field H s . Since LRO is a 3D phenomenon, we 
supply the 2D J\ — J[ — J2 model with a uniform inter- 
layer coupling Ji . The ordering transition is reflected by 
a change in the spin stiffness p s . For a finite lattice with 
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FIG. 4. (Color online) Temperature-vs-field phase diagram 
for Pb2Vs09. The solid and dashed lines show the temper- 
atures of the BEC (3D model) and BKT (2D model) transi- 
tions, respectively. Circles,— squares,— and triangles- 1 — rep- 
resent the experimental data. The shading denotes the region 
of the BEC phase. 



the fixed aspect ratio, the spin stiffness at the transition 
temperature Tbec scales as p s (Tbec) = L 2 ~ D , where 
D = 3 is the dimensionality of the systemj 27 i 28 Therefore, 
Lp s (TBEc) should not depend on L, and the transition 
temperature can be found as a crossing point of Lp s (T) 
curves calculated for different L. Since our system is 
strongly anisotropic, the aspect ratio of the finite lattices 
should be adjusted to achieve the proper scaling^ Fol- 
lowing Refs. [27] and l28l. we reduce the lattice size along 
Jl, and use 30 L x L x L/2 for 0.01 < J x /Ji < 0.05 as 
well as L x L x L/4 for J x /Ji = 0.005. 

Experimental information on Tbec of Pb2V30g is 
given by different techniques. The most accurate es- 
timates are available for (loH < 14 T from specific 
heat ££22 magnetic susceptibility^ NMR,^ 3 - and 
measurements. In higher fields, Tbec was probed 
by magnetic torque^ and susceptibility^ However, the 
present high-field measurements are subject to signifi- 
cant errors so that the estimates of the maximum Tbec 
range from 4.5 to 5.5 K in fields of 20 — 27 T (Fig. gj}. 
In our simulations, we first focused on the region below 
14 T and selected J±/Ji = 0.005 (about 0.15 K), which 
yields the best agreement with the experimental data be- 
low 14 T. The value of J± supports our DFT result of 
J± ~ 0.6 K. Note that a weak J± has nearly no effect 
on the magnetic susceptibility or magnetization curves, 
hence the interlayer coupling can be accurately estimated 
from Tbec only. 

The full temperature-vs-field phase diagram is shown 
in Fig. |31 Our simulations predict the maximum Tbec of 
4.4 K at p, H = 20 T in perfect agreement with the 
results of the torque measurements.— In contrast, the 
high-field susceptibility measurements 3 -* likely overesti- 
mate the transition temperatures. Specific heat or mag- 
netocaloric effect studies above 14 T should further chal- 
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FIG. 5. (Color online) Field-dependent temperature of the 
magnetic susceptibility maximum (Tm ax ) in PD2V3O9. Tri- 
angles represent the experimental data of Ref. |3lj. 



lenge our prediction. We also calculated the tempera- 
tures of susceptibility maxima (Fig. [5]) and found excel- 
lent agreement with the experimental data of Ref. [3ll 

To understand the origin of the BEC in Pb2V3 0g, we 
first calculate spin-spin correlations in zero field. We 
find that the system is close to the dimer limit, with the 
large correlation of —0.219 on the J\ bond and weaker 
correlations of —0.061 and 0.031 on the J[ and J2 (Jq) 
bonds, respectively. Magnetic field induces singlet-to- 
triplet flips on the spin dimers, whereas the interdimer 
couplings J[, J2, and J' 2 lead to the LRO. The reduc- 
tion in J 2 {J 2) suppresses the BEC (Fig. [HI top). The 
AFM coupling J[ has a similar, yet weaker effect (Fig. [51 
bottom). The difference should be related to the fact 
that the FM couplings J2 and J 2 join the dimers into a 
2D network, whereas the AFM coupling J[ connects the 
dimers into chains and does not lead to the LRO (the 
Ji — J[ model is ID). Based on these results, we identify 
the FM couplings J2 and J 2 as the main driving force of 
the BEC in Pb 2 V 3 0g. Note also that J 2 (J 2 ) provides 
an overall energy of — 0.5 Ji (two bonds per site) which is 
comparable to 0.65Ji arising from J[. This further jus- 
tifies the assignment of Pb 2 V30g to quasi-2D systems of 
coupled spin dimers. Despite leading to a reasonable fit 
of the magnetic susceptibility, the ID alternating-chain 
description is apparently incomplete, and does not cap- 
ture the essential microscopic physics of the system. 

The interlayer coupling J± is also important for the 
BEC. For example, J±/Ji = 0.01 results in the maximum 
Tbec of 4.8 K (compare to 4.4 K for J x /Ji = 0.005). 
The experimental data give a strong evidence for a non- 
zero Jx, and do not fit to a purely 2D model. In two 
dimensions, the LRO in a Heisenberg system is possible 
only at zero temperature, but the magnetic field induces 
an XY anisotropy giving rise to the finite-temperature 
Berezinsky-Kosterlitz-Thouless (BKT) transition. Using 
the scaling procedure from Ref. H3, we estimated Tbkt 
for the 2D J\ — J[ — J2 model.— The field dependence 
of Tbkt (dashed line in Fig. 2]) largely resembles that of 
Tbec j but the BKT temperatures clearly underestimate 
the actual transition temperatures in Pb2VaOg. A simi- 
lar effect has been reported for other quasi-2D systems: 
spin dimers on a square lattice^ and a simple square 
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FIG. 6. (Color online) The maximum BEC temperature de- 
pending on the FM interdimer couplings J2 and J2 (top) and 
the AFM interdimer couplings J[ (bottom). In the upper 
panel, open and filled circles denote two opposite regimes with 
J2 = J2 and J2 = 0, respectively. The interlayer coupling is 
Jx/Ji = 0.01. 

latticed Note that the frustration of interlayer couplings 
could lead to a peculiar high-field behavior in the vicin- 
ity of H c (Ref. 36). However, our DFT results suggest 
unfrustrated interlayer couplings in PD2V3O9. 

Finally, we explore the difference between J 2 and J 2 . 
To study its effect, we set J2/J1 = —0.5 and J 2 = 0, thus 
transferring the full energy of the FM exchange to the J 2 
bond while keeping the same J2/J1 = —0.25 [here, J 2 = 
(J2 + J 2 )/2]. No changes in the magnetic susceptibility 
or magnetization curves are found. By contrast, Tbec 
is slightly reduced (compare open and filled circles in 
Fig. [6]). The decrease in Tbec can be understood as the 
reduction in the coordination number (two for J2 = J 2 
and one for J' 2 — 0) and the subsequent enhancement 
of quantum fluctuations. Otherwise, the distribution of 
exchange energies between different bonds of the lattice 
has only a minor effect on the magnetic behavior as long 
as the quasi- 2D nature of the system is preserved. This 
situation is typical and can be compared to the effect 
of spatial anisotropy on a frustrated square latticed or 
honeycomb latticed 

The proposed magnetic model of Pb2V3 0g is sup- 
ported by a direct comparison to the experimental data, 
and discloses the microscopic physics of this material. 
Although the alternating Ji — J[ chains form a backbone 
of the spin lattice, the system should be rather viewed 



as 2D, since the purely ID model overestimates the spin 
gap and does not account for comparable interdimer cor- 
relations along a' and b 1 . The 2D model of coupled spin 
dimers provides a realistic description of Pb2V3 0g by re- 
producing the temperature dependence of the magnetic 
susceptibility, the field dependence of the magnetization, 
and the field dependence of the susceptibility maximum. 
To achieve an accurate description of the field vs. tem- 
perature phase diagram, the weak interlayer coupling Jj_ 
is required. 

The quasi-2D nature of Pb2V30g should inhibit the 
observation of effects that are typical for ID systems. In 
particular, manifestations of the Luttinger liquid physics, 
predicted for an isolated alternating spin chain^ are 
likely concealed by the BEC transition. The distinct 
features of Pb2V30g are the simple spin lattice, the 
pronounced two-dimensionality, and the combination of 
FM and AFM interdimer couplings. Most of the known 
BEC materials reveal a complex arrangement of inter- 
dimer couplings (e.g., in TlCuCls)^ or inequivalent spin 
dimers (e.g., in BaCuSi20e)^ thus realistic modeling is 
a tough problem. In contrast, Pb2VsOg demonstrates 
the BEC physics on a relatively simple and well-defined 
quasi-2D spin lattice and on an experimentally accessi- 
ble field scale. We expect that the recent progress in 
crystal growth2±22, will stimulate extensive investigation 
of Pb 2 V30g in fields up to H s . The quasi-2D nature 
of the spin lattice raises the issue of the actual dimen- 
sionality, and calls for theoretical study of critical expo- 
nents. We also note that our model is presently restricted 
to isotropic exchange couplings. Since S^VsOg reveals 
Dzyaloshinsky-Moriya couplings, 9,40 a sizable exchange 
anisotropy is also possible in Pb2VsOg and should be 
further probed by single-crystal ESR experiments. 

In summary, we have studied the electronic struc- 
ture of Pb 2 V30g and presented the microscopic magnetic 
model for the BEC in this quantum magnet. The com- 
pound can be understood as a quasi-2D system of spin 
dimers with both antiferromagnetic and ferromagnetic 
interdimer couplings. Quantum Monte Carlo simulations 
of the magnetic susceptibility, high-field magnetization, 
and BEC transition temperatures are in excellent agree- 
ment with the experimental data, and yield accurate esti- 
mates of individual exchange couplings. The simple and 
well-characterized spin lattice of Pb2VsOg makes this 
compound a convenient model system for future studies. 

We acknowledge fruitful discussions with Oleg Janson, 
Franziska Weickert, and Nicolas Laflorencie. A.T. was 
funded by the Alexander von Humboldt foundation. 
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